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\ We study the dynamics of impurity bound electrons interacting with a bath of conduction band 

electrons in a semiconductor. Only the exchange interaction is considered. We derive master 
equations for the density matrices of single and two qubit systems under the usual Born and Markov 
approximations. The bath mediated RKKY interaction in the two qubit case arises naturally. It 
leads to an energy shift significant only when the ratio (-Rt) of the inter-qubit distance to the 
thermal deBroglie wavelength of the bath electrons is small. This bath mediated interaction also 
has a profound impact on the decoherence times; the effect decreases monotonically with Rt. 

I. INTRODUCTION 

i-Sh . The dynamics of the impurity bound electrons in semiconductors have been studied previously in the context 
of population relaxationi2i£. Experimental studies of electron spin relaxation were central to understanding various 
mechanisms of relaxation of nuclear and electron spins in 31 P atoms in Si. Since then, theoretical work has exhaustively 
considered the effects of various interactions among electrons and atoms on the relaxation of populations in different 
quantum states. Various researchers have considered the effects of spin-orbit coupling, hyperfine, exchange, Coulomb, 
C$ ■ and dipolar interactions. Current interest in the unusually long relaxation times predicted and observed in these 
systems stems from possible applications in solid-state quantum computing. The seminal proposal by Kane^ to use 
31 P atoms embedded in Si crystal as qubits identifies one strategy for the realization of quantum computer hardware 
in the solid state. 

Recent investigations of decoherence mechanisms in these schemes assume the presence of a strong magnetic field 
and low temperatures, typically, a few mK; this prevents spontaneous spin flips by breaking the degeneracy of bound 
electron spins^SiLS. The dominant decoherence mechanisms in this case include the hyperfine and dipolar interaction 
of electrons with nuclei, and dipolar-dipolar interaction among the electrons. Recent calculations-? suggest that at 
low temperatures decoherence is dominated by spin diffusion induced by hyperfine and dipolar interactions of qubits 
with the spin bath of Si atoms. In principle, and to a high degree in practice, this mechanism can be eliminated by 
■ using ultra pure 28 Si crystals. But it would be desirable to operate solid state quantum devices at low magnetic fields 
and higher temperatures, and the price for this is the importance of additional decoherence channels in the system. 
CN| ' One of these is due to the interaction of qubits with a bath of conduction electrons. Studies of decoherence due to this 
channel have appeared in literature only recentlySiifl. The problem investigated by Kim et al& concerns the exchange 
interaction of a qubit with a spin polarized one dimensional stream of electrons. However, it only addresses the spin 
flip rates of the qubit, and uses a bath that is physically different from an unpolarized gas of conduction electrons 
typically found in a semiconductor. The latter type of bath is used by Rikitake et alm& who study the effects of 
decoherence on the RKKY interaction of two qubits in a bath consisting of a non-interacting degenerate electron gas. 

In this paper we consider a simple model, similar to that of Rikitake et alm&, in which the spins of electrons bound 
to donor atoms act as qubits and scatter the conduction electrons via the exchange interaction. For reasonably low 
donor densities, we show that the conduction electrons form a Boltzmann gas at all temperatures. Therefore, in a 
Kane type model only a classical distribution need be considered. The temperature is considered high enough that 
the ratio of bound to free electron density ensures that interactions among qubits are negligible compared to their 
exchange interaction with conduction electrons. For P atoms in Si, these assumptions are satisfied, in the absence 
of magnetic fields, for donor densities of the order of 10 16 cm -3 or lower. The nature of the bath has important 
consequences for the temperature dependence of decoherence. Thus the results in the present paper are in contrast 
to those in of Rikitake et at, despite similar master equations obtained in both. Furthermore, the present paper 
addresses the Kane model more concretely and makes stronger connection between the parameters of the equations 
derived and fundamental properties of semiconductors. 

In the following, we first present a full master equation for the density matrix of a single qubit and obtain an intuitive 
analytical result for the decoherence and relaxation times. The result is similar in form to the phenomenological result 
obtained by Pines et alSk for the relaxation rates under conditions similar to those considered here. We then derive 
the master equation for a system of two mutually non-interacting qubits, and study the decoherence due to their 
interaction with the conduction electron bath. We plan to include the effects of interactions among qubits in a future 
paper. 
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II. MODEL AND FORMAL EQUATIONS 

We consider a silicon lattice at nonzero temperatures doped with a density nr> of phosphorus atoms. Each P atom 
donates an electron, which either becomes a conduction electron or is captured by another ionized P atom forming an 
"atom" with hydrogen-like properties. The captured electrons are usually in s-states with a "Bohr radius" of about 
25 angstroms and a binding energy of about 0.044 eV&ii. The conduction electrons form a gas of approximately free 
particles with an effective mass of 0.2m e , where m e is the bare mass of the electron. The density of the gas builds 
up (from zero at T = 0) as temperature rises and more donors are ionized. A simple statistical mechanics analysis 
shows that the expected number density of bound electrons at temperature T is rib = + \ Z ~ X ex P( — Eb/ksT)]^ 1 7 
where z is fugacity of the total system comprised of bound and unbound electrons, Ef> is the energy needed to excite a 
bound electron into the conduction band, and kg is Boltzmann's constant. The conduction electron gas has a number 
density of 

n c (T) = 2X T 3 F 3/2 (z), (1) 

where F 3 / 2 ( z ) = dxypx{z~ 1 e c + l) -1 is the Fermi-Dirac function, and At = h(2ir/mkBT) 1 ^ 2 is the thermal 
deBroglie wavelength. Forcing n c + = hd , we find 

F V2 (z)[2z + e- E ^ T ] = ln D X 3 T e- E ^ T . (2) 

For the parameters chosen, the product X 3 r e~ Eb ^ kBT < 10~ 24 cm 3 for temperatures below 300 K, which implies that 
for 7i D ss 10 16 cm -3 , we can take z <C 1 and F 3 / 2 (z) ~ z ■ Consequently, the distribution of conduction electron gas 
remains Boltzmann down to T — 0. 

We simplify the picture by assuming the qubits to be in s-states, with only the spin acting as their degree of freedom, 
and suppose there are no external fields breaking the degeneracy of the spin states. The conduction electrons collide 
occasionally with the qubits elastically; we assume they do not excite them into higher states. However, the spins 
of the two may become entangled or even exchanged. Since the conduction electron moves throughout the crystal 
interacting with many electrons and atoms before the qubit scatters another conduction electron, it loses its coherence 
much faster than the qubits, and may be thought to be in an incoherent superposition of momentum eigenstates over 
the timescale of interest. Thus it is also independent of other bath electrons, while being on the same footing as them. 
This allows us to study only the case of interaction between a qubit and one mixed state conduction electron, and 
multiply by the number of the latter in the end result. Thus we can take our Hamiltonian as 



H = f- + V, 
2m 



(3) 



where p is the momentum of the conduction electron and m is its effective mass. The operator V is the interaction 
Hamiltonian acting on the electron and the qubit (s). We are concerned only with the exchange interaction, and study 
two cases. We first consider a single qubit at the origin for which the interaction takes the form 



V = Jr^(r)S-s, 



(4) 



where J is the exchange coefficient and ro is an effective Bohr radius characterizing the size of the qubit. The 
spin operators S and s act on the qubit and the conduction electrons respectively. We also consider two mutually 
non-interacting qubits for which the interaction is 



V = Jrt 



5(r - 1r)Si • s + 5(r + ^R)S 2 



(5) 



where Si and S 2 are the qubit spin operators, and the two qubits are placed symmetrically at ±R/2. To treat more 
than two qubits, similar terms would be added to V, with delta functions centered at the respective qubit locations. 
Despite the absence of direct exchange interaction between them, the two qubits can still exchange spins via indirect 
exchange interaction. Physically this occurs when the bath remains correlated long enough for the conduction electron 
to link two qubits; the indirect exchange coupling between qubits can arise even when they are too far away to have 
significant direct exchange interaction. This coupling is significant only when the inter-qubit distance is much shorter 
than the coherence length of the bath, which is approximately the thermal deBroglie wavelength At- The results 
found here depend naturally on these two important length scales. 

In addition, there are also direct interactions between the qubits. The important qubit-qubit interactions involve 
the exchange interaction between bound electrons, the hyperfine, dipolar, and spin-orbit coupling of these electrons 
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to the bath of Si nuclei. For low doping densities the first of these can be ignored as a starting point. For a donor 
density of ri£, = 10 16 cto~ 3 , and T = 100A", the inter-donor distance is R « 50 nm, whereas the mean radius of 
the electron orbits is tq w 2.5 nm. Thus the exchange energy is small, as we expect little overlap between the 
qubit wave functions. The second and third have been studied by various authors in the context of both relaxation 
and decoherence ratesi£*2i£*Li2iI£. The hyperfine and dipolar terms can be made arbitrarily small by purifying the 
Si samples to contain only the 28 Si isotope. For natural Si, which contains 95.33% 28 Si, it is estimated that the 
hyperfine and dipolar couplings combined are of the order of 10 eV or less 8 -. This is minute compared to the 
exchange interaction on the order of meV that we consider. Spin-orbit coupling is likewise small, and we neglect these 
additional effects in this preliminary investigation. 

We point out that the model has important differences to the one used by Chang et al. in their general study of 
dissipative dynamics of a two-state system^. They consider a biased qubit that is coupled to the bath via S z only, 
and the coupling induces no spin-flips in either the bath or the qubit. Furthermore, spin flips are introduced by a 
tunneling parameter that is independent of the bath state. This is clearly not the case in (|4I5|I where the isotropic 
coupling of system and bath causes joint spin flips in the two subsystems. A more general case of Brownian motion 
in a fermionic environment has also been studied in several papers by Cheni£*i&i£, who mainly focused on mapping 
between fermionic and bosonic environments. None of these studies discusses decoherence directly, and it is highly 
nontrivial to extend the results of these papers to arrive at those in this work. 

We first develop a general equation for the density matrix p(t) of a system coupled to a bath and the full system 
evolving via the Hamiltonian in ©. The following notation is used. We label the conduction electron states by 
| op) = | a) ® |p), where |a) and |p) are the eigenstates of s z and the momentum operator of the conduction electron 
respectively. The kinetic energy of the bath electron is denoted by huj p = p 2 /2m, and the normalized system states 
are labelled using Roman letters. For a general interaction V, we write the interaction Hamiltonian in the interaction 
picture as 

Hi(t) = E l*> 01 ®EE | ap ) ( ap | VtJ \a'p') (a'p'\ (6) 

ij a, a' p,p' 

Here we defined bath operators Vy = (* | V | j) which depend parametrically on system states. It is shown in the 
Appendix that for an unpolarized bath under the Born approximation, the density matrix p of the system is given by 

p = - iuirkky ^ T j j a [rjjj r)ki, p] + 7 A i j u {^Vki PVij - Vij Vkl P ~ PVij Vki), (7) 

ijkl ijkl 

where the operators r\ij — \i) (j are system operators. The tensors T and A, and the constants 7 and ojrkky, result 
from specializing (|A7|) and (|A11(1 to the two the forms of V shown in @ and J5J) . It is evident that T is involved with 
a unitary evolution of p(t) and A with the decoherence processes. Thus we label the first sum as the unitary term and 
the second as the dissipative or non-unitary. The frequency lorkky is a shift resulting from the RKKY interaction 
between the two qubits. No such term arises in the single qubit problem. The constant 7 -1 is an interaction time 
constant of the system and the bath, (see appendix for derivation); 7 and ojrkky are defined as 



7 = 4(27r)-ln c (T)(Jr W 2 )V^, (8) 

V m 

urkky(B) = -(27r) 3 h- 1 (Jr 3 ) 2 x(H) (9) 

= in c {T){Jrlmh- 2 f^-e-< 2R '^ 2 . (10) 
7r mK 

In the definition @ of ujrkky, x( r ) is the f ree electron susceptibility, and the expression in (I10|) specializes this defi- 
nition to a Boltzmann distribution for the conduction electron gas. Apart from dimensionless factors, 7 is a product of 
the thermal flux of conduction electrons and an "effective area" ( Jr^mh -2 ) 2 over which the bound electron experiences 
this flux. On the other hand, the thermal flux in the definition of ojrkky is replaced by n c h(mR)~ l exp(— AirR 2 ^), 
which depends on the inter-qubit distance and the thermal deBroglie wavelength of bath electrons. Both these rates 
are given by a mean number of scattering events, where the scattering occurs via the exchange interaction. The 1/R 
divergence of lurkky is typical of the x( r ) m three dimensionsi&iS but the spatial damping factor in the case of 
Boltzmann distribution is Gaussian as opposed to exponential in the degenerate limibiS,. Furthermore, the increase 
in decoherence with n c and J is accompanied by a correspondingly faster unitary evolution, as the ratio j~ 1 ujrkky 
is independent of both these parameters. 
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III. APPLICATION TO TWO SPECIFIC CASES 
A. Dynamics of a single qubit 

This is the simplest application of the general expressions given above, and we use it to estimate the temperature 
dependence of relaxation and decoherence times. The tensors T and A in this case are equal and their tensor elements 
are found to be 

A m = \ (i| S* \j) (k\ S* \l) + \ (i\ S + \j) (k\ S- |0 + \ (i| S- |j) (fc| S + \l) , (11) 

where S ± = S x ± iS v . Substitution of this expression in yields the master equation for a single qubit explicitly 
in the Lindblad form. The unitary term becomes [S 2 ,p], and it vanishes because S 2 is proportional to identity. We 
define a superoperator £(.) such that C{A)p = 2ApA^ — A^Ap — pA^A, where A is a general operator, and obtain 

p = l{£(S x )p + £(S y ) P + £(S z )p} (12) 

= -l([s*, [s*, P ]] + [s* [sy, P ]} + [s*, [s z , P ]]). (i3) 

The latter equality follows due to the Hermiticity of the spin operators. Besides ensuring complete positivity and 
trace preservation, this form is also invariant under arbitrary rotations of the spin. Therefore we may consider any 
pure initial state to be a spin up state in the S z basis of a suitable coordinate frame. In the Bloch formalism, these 
states reside on the surface of a unit sphere and their loss of purity is described by the decay in the length of the 
vector representing them. From l|13l) we then find that the initial decay of purity is equal to dTr[p 2 ]/dt = —7 for all 
pure states. 

To study the evolution of a general state, we study the Bloch vector with components u = poi + pio, v = i(poi ~ P10)) 
and w = pn — poo- Here identifies the "spin down" state and 1 the "spin up" state. From the master equation H13|) 
it follows that 

u = — 7U, v = — jv, w = —jw. (14) 

From these equations we find that the longitudinal and transverse rates are T-j - = T 2 _1 = 7. Thus the relaxation 
rate is proportional to the collision rate of a qubit with a thermalized bath of particles. The two rates are equal because 
the system is unbiased, and both relaxation and decoherence arise only through elastic scattering. The conventional 
relation 2T2 < T\ holds true only when several distinct scattering processes determine T2 while only a subset of them 
is responsible for T\ . When no such distinction exists we can have T\ — T2 , as Bloch noted many years ago^S. 

In Si:P with physical parameters defined in Section[nJ we estimate J k, 6meV. Taking into account the temperature 
dependence of n c (T), shown in Figure QJ, we plot log(7 _1 ) as a function of T in Figure (J2J) (the logarithm base is 
10). It is evident from the figure that above about 70 K, 7 _1 is less than a microsecond, which means that conduction 
electrons present a significant decoherence mechanism in this regime. The effectiveness of this channel vanishes 
significantly at low temperatures due to the loss of conduction electrons (as discussed in Section [HJ; n c is less than 1 
percent of ri£, for T < 35K. The y/T dependence of the thermal flux, which is associated with a decrease in conduction 
electron velocity with lowering temperature, also contributes to the rapid decrease in 7 as temperature is decreased. 

B. Dynamics of two mutually non-interacting qubits 

We now consider a system of two mutually non-interacting qubits that relax and decohere via the exchange inter- 
action with the conduction electrons. Again, the general dynamical map J7J describes the evolution upon specializing 
the tensors A and T to (JSJ) using the defining equations i|A7() and IjAllI) . As expected a contribution of the form 111|) 
for a single qubit comes from each member of the system. We call this A' and write it as 

2 1 1 1 

KjH = E 2 (?l S « W S « 10 + 4 M S « IJ> < fc l S « 10 + 4 M S * IJ> < fc l S « l*> ' ( 15 ) 

q=l 

where q = 1,2 labels each qubit. However, there also exist "cross-coupling" terms which represent the process by 
which a conduction electron mediates a spin exchange between the two qubits. These are given by A" where 

Kki = \ (»1 si \J) (k\ si 10 + \ (i\ S+ \j) (k\ 5 2 " 10 + \ ii\ S{ \j) (k\ s+ 10 
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Figure 1: Ratio of conduction electron density to donor density as a function of temperature 
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Figure 2: Decoherence time for a single qubit as a function of temperature (logarithm base 10). 
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\ (i\ S-2 \j) (k\ Sf \l) + \ (i\ S+ \j) (k\ 5f \l) + \ (i\ S 2 ~ |i) (k\ S+ \l) . (16) 



The two tensors add to form the tensors that appear in (JjJ) : 



A ijkl = A^ + £(i* T )A£ fci , (17) 
Tijki = A-ijki, (18) 



t{Rr) = xe-^[R T VA^]d X = Z^&^ e -^*l (19 ) 
Jo 4i?T 

where Rt = R/Xt is the inter-qubit distance measured in units of the thermal deBroglie wavelength, and Erf(.) 
denotes the error function. The tensor T depends only on A" because the term corresponding to A' commutes with 
p(t), as can be verified from Q. The dimensionless integral £(Rt) represents the strength of indirect exchange 
coupling between the two qubits relative to their direct exchange interaction with the conduction electrons. The most 
dominant contribution to the integral comes from x < ir/4R^. Thus we find, as expected, that the strength of indirect 
exchange decreases rapidly as qubits move out of the coherence region of the conduction electron, and the plot of 
£(Rt) m Figure Q shows a monotonic decrease as Rt increases. We point out that if the Fermi distribution were 
applicable, £(R,T) would be given by the same formula as above, but with xe~ x replaced by the x(z~ 1 e x + 1) _1 , 
where z is the fugacity of the gas. At vanishingly small temperatures, the Fermi wavelength would then take the role 
of the thermal deBroglie wavelength in setting the length-scale of indirect exchange. But we stress that for the model 
presented in section [n] the Boltzmann distribution is relevant for all temperatures. 

Complete positivity and trace preservation is guaranteed in the two qubit case as well. Substitution of I|15I16|I in 
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Figure 3: f (Rt) as a function of Rt = R/At 




yields 



P = 



-\Heff, P ] + c ^F kP F} - {F}F k ,p}], 



where we have defined an effective Hamiltonian 

Heff — ^RKKY^l ■ S2, 



(20) 



(21) 



which causes unitary evolution due to system bath interaction. The dissipative part contains six operators 
(F±, F 2 , F 3l F 4 , F 5 , F 6 ) — (Sf , S± , 5f , S%, S 2 , S 2 ), and a 6 x 6 symmetric coefficient matrix 



a = 4 



a £a 
£a a 

2 
1 
1 



(22) 
(23) 



The matrix c is positive semidefinite, as can be verified from its non- negative eigenvalues for < £ < 1. This is 
sufficient to ensure complete positivity of H20|) . The equation is reduced to Lindblad form by diagonalizing c and 
obtaining an orthonormal set of eigenvectors. Thus 



■j[%iP] 

i 7 (l 



sf) + r(sj + ^) + r(s| + sf)]p 

5f) + £(5|-^) + £(5|-5D]p. 



(24) 



An equation of the same form has been derived by Rikitake et al. but with the parameters specific to a degenerate 
bathiS. Let us now compare this map with the single qubit map (|12|l by ignoring the unitary term. For the case 
£ = 0, c is already diagonal, and the spin operators of each qubit form the set of Lindblad operators in the dynamical 
map. The map then consists of a sum of maps l|12|) for each qubit, which implies that the reduced dynamics of each 
qubit is independent of the other. Physically, at £ = 0, the coherence length of the bath is much smaller than the 
inter-qubit distance, and therefore the qubits scatter bath electrons independently of each other. 

In the presence of indirect coupling, the Lindblad operators are not the spin operators of the qubits but their sum 
S = S2 + Si, and difference A = S2 — Si. For significant values of £, the coherence length of the bath covers the 
inter-qubit distance. When scattering the bath electrons, it is then reasonable to expect that the two qubits behave 
as a single entity with total spin S. In fact, in the extreme limit of full coherence over the region containing the 
qubits, £ = 1, and the map has exactly the same form as that of the single qubit map. The difference operator, A, 
accounts for the deviation from perfect coherence of conduction electrons at the two sites, and allows independent 
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evolution to take place; the magnitude of this effect is of order 1 — £. This suggests that the singlet state with a total 
spin zero should be free of dissipative dynamics whenever £ = 1, and this is found to be the case in the solution of 
(124(1 presented below. Note that this is true not just within the context of the Born approximation to the scattering 
amplitudes. Since £ — 1 strictly only for R = (it is meaningless to consider T — as n c (0) = 0), the interaction 
becomes V = Jd(r)s ■ S. The operator S has a nullspace with the singlet as its only member, and therefore, the 
singlet stops interacting with the bath when £ = 1, and remains pure indefinitely. More generally, the purity is clearly 
long lived whenever At 3> R, and £ w 1. Note also that no state other than the singlet can be in isolation because 
a nonzero total spin would always interact with the bath electrons. Finally, the dynamics of the singlet was briefly 
considered by Rikitake et at. They find similar behavior for the singlet, but where the thermal deBroglie wavelength 
is replaced by Fermi wavelength, as was pointed out earlier for the case of a degenerate gas. 

Let us now consider in detail a realistic case of intermediate values of £. In order to proceed, we analytically solve 
(|24|l . which is done most conveniently in the Bell-basis representation of p(t), where the basis vectors are enumerated 
in the order they are shown below: 

[ 100) + jll) |00) - |01) + |10) |01) - |10) 

The first (second) symbol corresponds to the first (second) qubit. In this basis the operators S and A take a 
particularly simple form as shown in Table (JIJ. It is evident from the table and l|24|l that the diagonal terms of 
the density matrix do not couple to the off-diagonal ones. This simplifies the calculation and yields the following 
population equations: 

j t (p22-pn) = -1{2 + 0(P22-Pn), (25) 
^(p 33 -pii) = -7(2 + 0(P33-A»n), (26) 
j Y, Pa = -27(1-0 E + (27) 

i=l,2,3 1=1,2,3 

j fP44 = -27(l-0p44+|(l-0- (28) 

The first two equations show that the populations within the triplet manifold equilibrate to a uniform distribution at 
a rate of 7(2 + £). The last two show that the population transfer between this and the singlet manifold occurs at 
a rate of 27(1 — £), confirming our observation that the singlet ceases to evolve at £ = 1. The off-diagonal elements 
couple only to their conjugates. Within the triplet manifold they obey 

Pl2 {t) = (K[pi 2 (0)]+»[p 12 (0)]e-^ 1+ ^)e-^, (29) 

Pl3 (t) = ($t[p 13 (0)}+i%[p 13 (0)}e-^ 1+ &)a-^, (30) 

P23 (t) = (^f^[p 2 3(0)] e -^ 1+ ^ t -(-^3[ /32 3(0)]) e ^^ (31) 

while the elements between the triplet and the singlet manifolds obey 

pu(t) = ^[p 14 (0)]|cos(^)-^ + 7 0^^}e^ (3 -« )t/2 

+»9[M0)] {cos(c't) - (iw - 7')^^} e -^-^l\ (32) 
P2,(t) = ^[ P24 ( )]|cos(^)--^-70^^}^ 7(3 - €)t/2 

+^[ P24 (0)] jcos(^) - (iw + 7')^^} e -7(3-C)t/ 2j (33) 

P3A(t) = ^[p 3 4(0)]|cOs(^)-^- 7 0^^}^ 7(3 - €)t/2 

+l 3[p 34 (0)] f coB(w'i) - (iw + 7')^^) e~^-^l 2 . (34) 
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Here we defined the difference of eigenvalues of H e ff in the triplet and singlet manifolds as Hut = {k\H e ft \k) — 
(4|iJ e //|4) , where k = 1,2 or 3 labels the triplet states. The renormalized frequency u>' = y 'w 2 — r y' 2 , where 
7' = 7(1 — £)/2, represents the oscillations caused by the unitary evolution resulting from the RKKY interaction, ft 
follows from Table |QJ that these oscillations are absent within the triplet manifold; they also disappear for uj < 7' in 
the above equations. We note that the off-diagonal elements always decay at least with a rate of 7, and each of these 
elements is uncoupled from all others. Hence dephasing between any pair of Bell states proceeds independently of the 
rest of the states. 

Equations I|25I28[) show that, for £ 7^ 1, the final state of the density matrix is the maximum entropy state p = hi. 
However, for £ = 1, the relaxation between the singlet state and the triplet manifold ceases, and the final state 
becomes p{™ al = /344(C)) and pf/" ai = |(1 — /944(C))) for i = 1, 2, 3. The relaxation rate of the singlet can become zero, 
but all other states attain a minimum relaxation rate of 27. 

Several other properties of the dynamical map (|24(l become evident when we consider the rate of decrease in purity 
of an initially pure state; the map ensures that the purity p{t) = Tr[p 2 {t)] is a monotonically decreasing function of 
time. A general equation for the rate of change of purity follows straightforwardly from l)24|) as 

J = -3 7 (1 - £)p " 2 7 £Tr[SV] + 7 Tr [(1 + £)(Sp) • (Sp) + (1 - £)(Ap) • (Ap)] . (35) 
For pure initial states p(0) — |V) (V|, the initial loss in purity occurs at the rate 

p(0) = -37(l-£)-27£(V;|5 2 |^+7(l + 0ll(^|S|^)|| 2 +7(l-0ll(^|A|^)|| 2 . (36) 

ft is easily verified that the separable states of the general form \a) \b) lose initial purity at the rate p(0) = j(— 3 + 
2||(a|Si|a)|| +2||(6|S 2 |6)|r), which is independent of £. Thus the size of £ has no effect on the initial decay in 
purity of unentangled states. Furthermore, since any state corresponds to "spin up" in some direction, we have 
(a\ Si I a) • (a| Si |a) = (b\ S2 \b) ■ (b\ S2 \b) — 1/4, and all separable states are on an equal footing with respect to the 
initial rate of loss in purity. 

We now show that this is in fact the minimum rate achievable within the triplet manifold. An arbitrary state in 
the triplet manifold has (V| S 2 \ip) — 2, and (tp\ A \ip) = 0, which when substituted in l)36|l yields the initial rate 
p(0) = -7(3 + £ - (1+ || (VI S |V) II ) with only one state dependent variable, ||(V| S |V)||- The maximum value 
of ||(V|S|V)|| occurs for a separable tp, which consists of parallel spins. Therefore, we see that within the triplet 
manifold, separable states are the most robust against loss in purity. 

On the other hand, the singlet has \p(0) \ — 37(1 — £), which becomes less than 27 for £ > 1/3. Therefore the singlet 
is more robust than the separable-states in this regime. The rate corresponding to any superposition of singlet and 
separable states is greater than the average of the rates of these states. Consequently, the set of robust states does not 
change in any continuous manner with £, and it instead changes from separable to singlet at £ = 1/3. For t > 0, the 
fidelity, f(t) = Tr[p(i)p(Q)], offers a characterization for predictability in which f(t) as 1 identifies a highly predicable 
state and f(t) ~ a poorly predicable one. In a study of fidelity, not reported here, we found that the set of most 
predicable states makes a transition from separable to singlet state for 1/3 < £ < 1/2, depending on the time elapsed. 

We end this section by discussing the time-dependent behavior of p(t) for a few simple cases. Since for our Hilbert 
space of dimension 4, p(t) decreases monotonically to the value 1/4 for £ ^ 1, it is useful to define an "instantaneous 
rate" for purity 



a(t) 



P 



p-0.25' 



in terms of which 



p(t) = I + ^exp(^-^ a(t')dt' 



The function a(t) highlights differences in p(t), and it is especially useful when different states lead to similar behavior 
for p(t). We first consider the separable states of the form 

|V) = |O)(cos(0)|O)+sin(0)|l)) 



Table I: Bell basis representation of the spin operators S and A. 



j -» 


X 


y 


z 




|1> (3| + |3) <1| 


-i\2) <3|+z|3) (2| 


|1> (2| + |2>(1| 




-|2) (4|-|4) <2j 


<4|+z|4)(l| 


|3> (4| + |4> (3| 
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Figure 4: The "instantaneous rate" a(t) for (a) £ = and (b) £ = 0.7 




(a) (b) 



and plot the ^-parameterized rate a$(t) as a function of <j> and time t for different values of £. Figure shows that 
all separable states give rise to exactly the same rate at all times for £ = 0. Here we set to = to describe the effects 
of non-unitary dynamics only. The effect of non-zero exchange coupling becomes evident in Figure flip) where we set 
£ = 0.7. Here the function a$(t) decays much more slowly with time for states <j> « 7r/4 than for <j> £ {0, 7r/2}. Thus 
while £ plays no role initially for the separable states, the more stable states are those with both qubits prepared 
parallel or anti-parallel to each other (see discussion after (13611 ). Similarly, the most vulnerable of separable states are 
those in which the qubits are eigenstates of spin operators corresponding to orthogonal Cartesian directions. 
The purity of the four Bell states has the following time dependence: 

Pi(t) = i + ie-^-Ot + He-** 2 *)*, i= 1,2,3 
Mt) = \ + \e-^-&. 

Thus the three states with total spin S ^ lose purity for all values of £, and do so at two different rates for £ ^ 1. 
Initially, the decay is dominated by the rate 27(2 + £), while at times much longer than the inverse of this rate, the 
decay approaches the slower rate 47(1 — £). The corresponding rate a(t) decreases from 47(1 + £/3) to 47(1 — £) as 
time increases. In contrast, the rate for the singlet is independent of time. 

IV. CONCLUSION 

In summary, we have derived and studied the master equations for impurity bound electrons (qubits) scattering a 
bath of conduction electrons in a semiconductor. We show that the distribution of bath electrons remains Boltzmann 
at all temperatures, due to the temperature dependent filling of the conduction band. Thus our analysis, based 
on this result, is in contrast with other studies on decoherence of qubits in an electron gas obeying a Fermi-Dirac 
distribution*!*^. The master equations are obtained in the Lindblad form for a single qubit and a system of two 
mutually non-interacting qubits. In the former, the Lindblad operators are found to be the spin operators for the 
qubit. In the latter, these are replaced by the sum and difference of the spin operators of each member of the system. 
The Bloch equations derived for a single qubit show that decoherence occurs at the same rate as relaxation in the S z 
basis. This departs from the conventional inequality 2T2 < T\ because the two levels in the system are degenerate, 
and therefore relaxation and dephasing processes are both elastic. The inverse relaxation times are equal to 7, which 
is proportional to the product of thermal flux of the bath electrons and an effective cross section that depends on the 
exchange coefficient J and the Bohr radius of the qubit. Calculations show that 7 _1 is on the order of seconds for 
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temperatures below 10 K, and decreases rapidly to below a microsecond above 70 K. Thus the conduction electron 
bath is dominant in causing decoherence compared to other sources^^LS for temperatures above 70K. 

The same parameter 7 also sets the timescale of decoherence and relaxation in a two qubit system. However an 
additional parameter, £(Rt), representing the indirect exchange coupling of the two members of the system, affects 
the rates profoundly. The function £(Rt) decreases monotonically from one to zero with Rt, the ratio of inter-qubit 
distance to thermal deBroglie wavelength of conduction electrons. The dissipative part of the dynamical map is 
dominated by the total spin operator when £ w 1. As a result the singlet state, with zero total spin, becomes the 
most robust state in this limit. For small £, however, the separable states in which both members of the system are in 
an eigenstate of the same Cartesian component of the spin operator exhibit the slowest rate of loss in purity. Pairs of 
Bell states are found to dephase independently of each other, and their dephasing rate never exceeds the population 
transfer rate. The unitary RKKY interaction arises naturally between the two qubits in our master equation. The 
frequency shift associated with this interaction is found to be proportional to the free electron susceptibility of the 
bath, in agreement with past studies of RKKY interaction between two spins mediated by a gas of free electrons 18 - 19 . 
While these studies found the interaction to decay exponentially as a function of inter-qubit distance for a Fermi-Dirac 
distribution, we find a Gaussian decay for a Boltzmann distributed electron gas. The results in the two qubit system 
can be understood qualitatively in terms of the coherence length of the bath electrons and the initial entanglement 
between the qubits. Electrons with large thermal deBroglie wavelength tend to scatter as if the two qubits were acting 
as a single entity, while those with small wavelength scatter off each qubit independently of the other. Similarly, qubits 
prepared in pure separable states lose purity independently of £, while the sensitivity to £ is much greater for entangled 
initial states. 

The generalization to include the effects of an external magnetic field is straightforward and will restore the in- 
equality 2T2 < T\ in addition to the precession of qubits in the field. However, interactions among the qubits demand 
a more involved calculation, because unless they are much larger or much smaller than the system-bath interaction, 
they render the secular approximation invalid. This approximation is central to most derivations of a coarse-grained, 
Markovian master equation. 
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Appendix A: DERIVATION OF MASTER EQUATION 

Here we derive the master equation {7J of section 2. The density operator, denoted W(t), for the joint qubit-electron 
system evolves unitarily via the transformation W(t) — U(t)W(0)U'(t), where the unitary operator U = 1 + iT and 
T satisfies the equation 

T = -±J HrfW-^J HzWTtfW. (Al) 

By unitarity of U(t) it follows that i(T - T+) = -T^T, and the evolution of W(t) can then be written in terms of T 
as 

W(t)-W(0) = ^[T + T\ W(0)] + TW(0)T^ — ^T^TW(Q) — -W(Q)T^T. (A2) 

We consider only the product initial state, W(0) = po <8> Rb, where Rb is the thermal density matrix of the bath. 
We first show that under the second order Born approximation and appropriate coarse-graining of t, (|A2fl yields the 
following equation for the system density matrix: 

p(t) - po = -i tNjjki [rnjilki , Po] + tM ijM (2r)ijp r) k i - rjkiVijPo - PoVkiVij)- (A3) 

ijkl ijkl 

The tensors Nijki and Mijki are independent of time, and 77^ are system operators defined by rjij = \i) 

We derive each sum in l|A3(l from the corresponding term in ljA2fl . correct to second order. Since the expansion of 
T starts at the first order in if/, the dissipative term becomes second order automatically. Application of the first 



11 



order expansion of T then yields the second sum in i|A3|l with the tensor Mijki given by the expression: 

tM m = ^_ fdt 1 /' t dt''^^e i ^-^'» t ''- t '^Qp;i|y|aV;j)(Qy;fe|y|ap;/)^p, (A4) 

a, a' p,p' 

where A/" p is the number of particles in the momentum state p. We now evaluate the integrals and sums in the limit 
of a large crystal. Let Q represent the crystal volume, and within this volume, let n(p) be the density of conduction 
electrons in the phase space. Then JV P = SlApn(p), where Ap = (2ttH) 3 /fl is the volume in p-space associated with a 
momentum state. Introducing the phase space volume for p' in a similar manner and letting f2 be large we find that 



E-^pE / d Pttn(p) J dp' 



n 



Next we introduce the rescaled operators Vij — £1{2-kK) 3 (i\ V \j) and substitute them in (|A4|) . Making a change 
of variables in t' , t" we find that 

ft r+L(t;r') r 

tM m = 47r 3 frV / dr' / dr / dpdp'e^-^ T n{p){ap\V lJ \a'p') {a'p'\V k i\ap) , 

aa , Jo J-Ht-T') J 

where L(t;r') = 2r' for t' < t/2 and L(t;r') = 2(r — r') for r < i/2. The integral over r is essentially a Fourier 
transform of a p~ dependent function and it is expected to drop quickly for thermal distributions at sufficiently high 
temperatures^. Therefore, we switch to a coarse-grained timescale which captures the evolution of the qubit but not 
the behavior of the bath correlations. In this limit, Lit; r') may be considered infinite for essentially all r' without 
introducing error in the integral over r. The result is 



tMijki 



aa' 



The factor t appears from the integral over r' due to the absence of system dynamics; in the presence of system 
dynamics, r' is coarse-grained further to be insensitive to the energy difference between the system levels. We now 
integrate over p' making use of the formula n(p)o?p = 4irp 2 n(p)dp(du/4ir), where n{p) is the occupancy of energy 
state p 2 /2m divided by phase space cell volume, and du is an element of the solid angle centered about u. After 
substituting these definitions we find that 



M ijU = 2{2i,f{mhf j°° dp^p 2 n{p) V / (apu\ % \a'pu') (a'pu'\ V u \apn) . (A6) 

Jo m ,J 47T 47T 

aa' 

In a similar treatment^ it was shown that a more accurate calculation yields a result obtained by replacing the first 
order scattering amplitudes in this formula with their exact counterparts. We expect the same to hold here, but 
since the exchange interaction is often introduced with parameters assumed appropriate for only a lowest order Born 
calculation, we do not pursue this issue here. We finalize this formula by substituting the Boltzmann distribution for 

n(p), 

n{p) = n c (2nmk B T)-h- p ^ 2mkBT , 

where we remind the reader that n c is the total number of conduction electrons present at temperature T. It is con- 
venient to introduce the dimensionless variable x — (2mkBT)~ 1 p 2 in terms of which (v — 4(27r) 9 / 2 n c (mS) 2 i/fc^T/m) 

M m = v J dxxe~ x J ^^E^' pfi '^J ' (a' ,pu'\V a \a,pu) . (A7) 

aa' 

Our next task is to expand iTrg (T+T')/2 and obtain a time- independent expression for Nijki ■ As the bath distribution 
does not depend on spin, the first order term vanishes by the zero trace property of Pauli matrices. The second order 
term then yields, 

~Tr B [T + T\W ] = ^IJ^dt'J^ dTTr fl [[fl J (t'),flj(*' ' - t)], po ® R B ]. (A8) 



12 



The commutator [Hi(t'), Hj(t' — r)] takes the following form where the summation is done over all p and a: 

[Hi(t'), Hi(t' - t)] = rji^ki J2 e l K"-"p)*' f e *(<v-«»)r _ e -i(u p ,-uj p „)A 
x |ap) (ap| Va \a'p') (a'p'\ V a |a"p") <a"p"l . 
We convert the summations over momentum to integrals and substitute the result in (|A8|) . Taking the trace we find 

l -Tr B [T + T\Wv] = [mjVki,Po]t(~^ 3 h)J2 J dpdp' 

a, a' 

poo 

/ dT2ism[(oj p , - uj p )t} (ap| % \a!p') (a'p'\ 7 W |ap) n(p). (A9) 
Jo 

The integral over r yields 1iP{lo v — ujp')" 1 , where V denotes the principal value. Comparing l|A9j) with (|A3|> . we get 

dpdp' n(p) V (ap| ^ |a'p') (a'p'| ? w |ap)]. (A10) 

ay — ui p — ' 



a, a' 



Substitution of the Boltzmann distribution and re-expression in terms of the dimensionless variable defined above 
yields the following expression: 

N ijkl = - [ dxdy^e-*— f V (q P | % \olp') (a'p'\ V kl \ap) . (All) 

IT J V — X J 47T 4-7T 

a, a 

Having shown the validity of l|A3(l . we can obtain a coarse-grained differential equation for pit) by iterating this 
equation after replacing p(t) — p(0) by pit + St) — p(t). Thus 

^ = Nijki[VijVkh p] + ^ Mjjki{2rikipr)ij ~ VijVklP - PVijVkl)' ( Al2 ) 

When the expressions and are substituted for V, the tensor M becomes 7A, where 7 is a constant defined 
in (JSJ) and A is given by i|ll|) for the single qubit and by (|17|) for the two qubit system. The commutator associated 
with the tensor N vanishes for the single qubit. In the two qubit system, the tensor N = ujrkky^, where ojrkky is 
defined in 10, while T is given by (|18f) . 

We now outline the calculation to obtain the RKKY splitting in terms of the susceptibility of the bath. We first 
observe that (|A10|I can be written as 

N m = 4n 3 h f dpdp'2^ ( n(p) V (ap| V l3 \a'p') (a'p'\ V kl \ap)}, 

J \ijJp' — Up + ley — ' 

where e — > at the end of the calculation. Substituting @ in the above expression, keeping only the cross-terms, and 
changing the variables of integration to k = p/h, wc find that 

N I1U = r, Ju .A-^f/A^„ w ,<.-» !8 (_^) +[ . t . 

Here k, q are wave- vectors, n(k)dk now represents the density of electrons with wave-vector within dk of k, and 
^ijki represents the summation over spin indices a, a'. We now manipulate the sum by first writing it as a sum 
of two identical copies of itself and then interchanging k, k' in one of the integrals. Then doing the transformation 
k, k' — > — k, — k' in that integral, we find that 



N ijk i = T ijkl 4ir 3 h- 2 (Jr 3 ) 2 



J {2-k) 6 [2Ti) i \aj k , ~ u 



- ui k + te 



^^ e .<^.*„ (k)s (_i_ 1 + ,, 
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We do not change the sign of ie in the second integral since the real part is unaffected by it. An expression equivalent 
to the one above is 

/ >j.r f dk [n(k+±q)-rc(k-§q )] 



N ijkl = -Y ijkl ^h- l {Jr 6 Q Y / ^e'"'" / - L v ^' v . 2 " /J +cc. 

J (2ir) 3 J (27r) 3 fi 2 k • q/m + 

The integral over k in the limit e = is the static Lindhard function 23 , which defines the Fourier transform, %(q), of 
the static susceptibility ^(r) 19 i 24 . Since x(r) = %(— r), 



JVy H = -(2 7 r) 3 ^- 1 (Jr3) 2 X (R)r 
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